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^^ , Abstract. We review VLBI observations of supernovae over the last quarter 

century and discuss the prospect of imaging future supernovae with space VLBI 
in the context of VSOP-2. From thousands of discovered supernovae, most of 
them at cosmological distances, ^50 have been detected at radio wavelengths, 
most of them in relatively nearby galaxies. All of the radio supernovae are Type 
II or Ib/c, which originate from the explosion of massive progenitor stars. Of 
these, 12 were observed with VLBI and four of them, SN 1979C, SN 1986J, SN 
1993J, and SN 1987A, could be imaged in detail, the former three with VLBI. 
c"j In addition, supernovae or young supernova remnants were discovered at radio 

Q_|. wavelengths in highly dust-obscured galaxies, such as M82, Arp 299, and Arp 

220, and some of them could also be imaged in detail. Four of the supernovae 
so far observed were sufficiently bright to be detectable with VSOP-2. With 
VSOP-2 the expansion of supernovae can be monitored and investiated with 
~* ' unsurpassed angular resolution, starting as early as the time of the supernova's 

transition from its opaque to transparent stage. Such studies can reveal, in a 
movie, the aftermath of a supernova explosion shortly after shock breakout. 



C ■ 1. Introduction 

q 

t — ■ A supernova (SN), the explosion of a star, is one of the most energetic single 

events in the universe. Thousands of optical SNe are now known but most 
of them are at cosmological distances. Only ~50 have been detected at radio 
wavelengths and 12 observed with VLBI, the most distant at almost 100 Mpc. In 
addition, several SNe and supernova remnants (SNRs) were discovered at radio 
wavelengths in dust-obscured galaxies, such as M82, Arp 299, and Arp 220, and 
observed with VLBI. All of the SNe observed with VLBI are thought to have 
resulted from the core collapse of a massive progenitor and emit synchrotron 
radiation generated from the electrons accelerated in the region where the shock 
front interacts with the circumstellar medium (CSM) left over from the wind 
of the mass-losing progenitor. In addition, some of the radiation may come 
from the environment of the stellar corpse, a neutron star or a black hole left 
over from the explosion. During the early stage of the expansion of a radio 
SN, the radio lightcurve rises quickly because of the decreasing optical depth 
due to synchrotron self-absorption or to thermal absorption in the ionized CSM 
along the line of sight. It reaches its peak, first at high frequencies later at 
lower frequencies, when the optical depth decreases below unity, and it declines 
there after during the optically thin stage of the SN while the SN is expanding 
fe.g.. lChevalienll982l b 

For a typical shock front expansion velocity of 10,000 km s , the angular 
size of a SN at a distance of 4 Mpc expands at a rate of 1 mas yr _1 . With global 
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VLBI and a wavelength of 4 cm, an angular resolution of ~0.6 mas can be ob- 
tained, allowing detailed investigations of the expanding SN and the making of a 
movie. Global VLBI with VSOP-2 promises to increase the angular resolution at 
this wavelength to 0.2 mas and to a correspondingly higher resolution at shorter 
wavelengths, enabling us to witness the aftermath of the explosion shortly after 
shock breakout. 



2. Why supernova VLBI? 

The high angular resolution of VLBI affords us the following possibilities: 

1. The determination of the morphology. A shell morphology indicates the 
interaction of the ejecta with the CSM. A centrally-condensed morphology can 
indicate the presence of a central source, such as a pulsar wind nebula or a black 
hole accretion disk system. A disk-like morphology suggests an optically thick 
object, as is expected in the early stages of evolution. 

2. The measurement of the angular size, 9, which can give us the average 
expansion velocity. It helps unravel whether synchrotron self-absorption played 
a dominant role in the early evolution of the SN, and could reveal any possible 
presence of a relativistically expanding afterglow of a gamma ray burst (GRB). 

3. The measurement of the expansion as a function of time, t, since shock 
breakout, parametrized by 6(t) oc t m ^' , which can give us the SN's age and 
deceleration parameter, m(t). 

4. Monitoring the expansion of the SN, which provides information on the 
density of the ejecta and the CSM, and on the magnetic field, as a function of 
radius. Obtaining the density profile of the CSM is like a "time machine" that 
reveals the last thousands of years of the mass-loss history before the star died. 

5. A geometric determination of the distance to the SN and its host galaxy using 
the ESM (expanding shock front method). 

6. The making of a movie showing the evolution of the SN. 
3. Supernovae observed with VLBI 



Supernova VLBI started with observations of SN 1979C in 1982 (JBartell [19851 ; 



Bartel et al.lll985l ) shortly after the new, sensitive, MKIII VLBI data acquisition 



system became available. Since then 11 more optically identified SNe have been 
observed with VLBI (Table [1]). Eight are Type II, whose progenitors are massive 
supergiants. Four are Type Ib/c, whose progenitors are believed to be the cores 
of massive stars which have lost their envelopes. A few Type Ib/c SNe were 
found to be associated with GRBs. The sources in M82, Arp 299, and Arp 
220 are not optically identified but are likely also SNe of Type II or Ib/c or 
SNRs. The distances range from 50 kpc for SN 1987A to almost 100 Mpc 
for SN 2003L. Their peak flux densities at 8 GHz reached values of 100 mJy 
and more, corresponding to very high brightness temperatures. Such SNe are 
examples for future targets of VLBI with VSOP-2. 



Supernova VLBI 3 

About 80% or more of all VLBI observations focussed on just three SNe: 
SN 1979C, SN 1986J, and SN 1993J. They were sufficiently bright and extended 
that detailed images could be obtained and, for the latter two, movies made. 

Table 1. Characteristics of supernovae observed with VLBI 



Name 


Type 


Host Galaxy 


Distance 


S 11 


< 2 t 


o 3 , 


Reference 4 








(Mpc) 


(mjy) 


(yr) 


(mas) 




SN 1979C 


II 


M100 


15 


7 


0.9 


0.30 


1 


SN 1980K 


11 


NGC 6946 


6 


2 


0.27 


<0.21 


2 


SN 1986J 


11 


NGC 891 


10 


100 


2.7 


1.6 


3 


SN 1987A 


II 


LMC 


0.05 


80 


0.005 


1.5 


4 


SN 1993J 


11 


M81 


4 


100 


0.27 


0.55 


5 


SN 1996cr 


11 


Circinus Gal. 


4 


150 


2.7 


2.8 


6 


SN 2001em 


Ib/c 


NGC 7112 


80 


2 


2 


0.11 


7 


SN 2001gd 


II 


NGC 5033 


13 


4 


0.5 


0.16 


8 


SN 2003L 


Ib/c 


NGC 3506 


92 


3 


0.2 


0.01 


9 


SN 2004et 


II 


NGC 6946 


6 


2 


0.13 


>0.14 


10 


SN 2007gr 


Ib/c 


NGC 1058 


10 


<1 


? 


? 


11 


SN 2008D 


Ib/c 


NGC 2770 


27 


3 


0.04 


0.01 


12 


SN/SNR 


1 


M82 


1 


1 


? 


7 


13 


SN/SNR 


? 


Arp 299 


-10 


1 


? 


? 


11 


SN/SNR 


? 


Arp 220 


77 


? 


? 


? 


15 



Notes: 

1 The approximate peak flux density of the radio lightcurve measured at 8 GHz or extrapolated from 
lower frequencies. 

2 The approximate time after shock breakout when the 8-GHz radio lightcurve peaked, i.e., when the 
optical depth decreased to below unity. If t pea k was not measured at 8 GHz, it is assumed to be 
0.7 X t pea fe at 5 GHz. For items 1 and 2, see references (last column) and/or K. Weiler's home page at 
http://rsd-www.nrl.navy.mil/7213/weiler/sne-home.html 

6 The approximate size of the SN at time £ pea fc extrapolated from the size and expansion given in the 
references. For SN 1987A, an expansion velocity of 35,000 kms" 1 is taken (Gacnslcr et al. 2007). For 
the Type Ib/c SNe, for which only upper limits of sizes and expansions were available, an expansion 

velocity of 10,000 kms -1 is assumed. 

4 Latest, or lates t comprehensive publications: 1) Bartcl & Bictcnholz (20C J3); 2) IBartell 1119851): 

3) Bictcnholz ct al. ( 2004^ 4 ) IJauncev et al.l (119881); 5) IBartel et al.1 1120001); 6) IBauer et alj 

(2008); 7) Bictcnholz & Bartcl (2007); 8) Perez- Torres ct al. (2005); 9) Sodcrbcrg ct al. (2005); 10) 
IMarti-Vidal et al.l (120071); 11 ) IParagi et alj 1120071) ; 12) Sodc rberg et alj 1120081 ); 13) IBeswick et alj 
(2006); 14) Neff et al. (2004); 15) Lonsda le et alj 12009) ■ Also, see these references or references in 
there for the characteristics in the table. 

In the remainder we will first focus on these three SNe. Then we will 
discuss the other optically identified SNe and the unidentified SNe/SNRs and 
lastly address prospects for space VLBI in the context of VSOP-2. 



4. Supernovae with detailed VLBI images 

4.1. SN 1993J 

SN 1993J in M81 is, after SN 1987A, the most intensely observed SN ever. It 
provides an exemplary case to witness, with increasing relative angular resolu- 
tion, the evoluti on of a SN over time . VLBI observations started ~ 30 d after 
shock breakout ( Marcaide et al.1 Tl994h 
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.1994 ) and showed that the 
SN expanded almost freely early on (jBartel et al.lll994l ). A first detailed image 
was obtained only seve ral months later after the SN had expanded sufficiently 
( Marcaide et alJll995al . see also, Bartel et al. 1995). Fig. [1] shows an example of 
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a sequence of images of the expanding shell of the SN ( Bietenholz et ahll2003l . 
see Marcaide et al. 1995b for an earlier sequence). 
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Figure 1. A sequence of 8.4-G Hz VLBI images of SN 1993J from t = 50 
to 3345 d after shock breakout (|Bietenholz et all I2003T ). The contours are 
at -1, 1, 2, 4, . . . , 32, 45, 64, and 90% of the peak brightness. For im- 
ages here and hereafter, the scale on the vertical axis is the same as that 
on the horizontal axis, the restoring beam is plotted in the lower left, and 
north is up and east to the left. (To download SN VLBI movies, please visit: 
http://www.yorku.ca/bartel.) 



The explosion occurred at a location pinpointed with 160 AU accuracy 
in the galactic reference frame of the host galaxy M81. Fro m there, the SN 
expanded isotropically to within 5.5% ( Bietenholz et al.ll200ll ). 

Initially, the SN expanded rapidly with a velocit y of ~ 20,000 k ms -1 . Then 



it was found that the expansion decelerated slightly (IMarcaide et al.l 11997) and 



that the deceleration changed over time, t, (Fig. 2, left panel; iBartel et al.ll2000l . 
2002) with the velocity decreasing to ~10,000 kmsj . These measurements are 



fairly consistent with a hydrodynamic model ( Mioduszewski et al.l l200ll ) and 
provide insight into the physical interpl ay between the ejec ta and the CSM, and 
the massdoss history of the progenitor ( Bartel et al.l 12002). 
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SN 1993J provided also the best example of geometrically determining a dis- 
tance using the expanding shock front method (ESM). The distance is derived to 
be 3.96±0.29 Mpc (Fig. [21 right panel) , which, when co mpared with Cepheid dis- 
tances, is somewha t larger than Freedman et al.l ([1994]) 's value of 3.63±0.34 Mpc 
but comparable to lHuterer et al.l ( 19951 ) 's value of 3.93±0.26 Mpc. 
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Figure 2. Left: The expansion of SN 1993J. The straight lines give fits 
which show the changing deceleration of the expansion with 9 oc £°- 919 ± - 019 
for the solid line and 9 oc £°- 781 ± 0009 f or the dashed line. Th e inset gives 
a zoo med-in version of the latest part of the expansion curve (|Bartel et al.l 
2002). Right: The velocity of the forward and reverse shocks (outer and inner 
radio shell) with their mean fit to the maximum Ha line velocities, giving th e 
distance to SN 1993J and M81 as indicated in the figure ([Bartel etal.ll2007| ). 



4.2. SN 1986J 

SN 1986J in NGC 891 is an intriguing extragalactic SN, the first one for which a 
compact component was found to emerge in the projected center of the expand- 
ing shell, ~20 years after the explosion. It was als o the first optically identified 
SN for which a detailed image could be obtained ( Bartel et al.lll99ll ). A series 
of images is shown in Fig.[3[ The new component has a spectrum which was 
inverted between 5 and 20 GHz when first discovered. The peak is slowly mov- 
ing to lower frequencies. The inverted part of the spectrum can be interpreted 
as being due to free-free absorption within the SN shell with the absorption 
weakening as the shell expands. The new component could be radio emission 
associated with accretion onto a black hole or the nebula formed around an en- 
ergetic young neutron star in the center of SN 1986J, which would be the first 
direct link of a black hole or a neutron star to a modern SN. Its spectral lumi- 
nosity between 14 and 43 GHz is ~ 200 times that of th e Crab Nebula, whic h 
is in the range expected for young pulsar wind nebulae ( Bandiera et al.|[l984l ). 
However, since the new component is in the center only in projection, it could 
also be caused by the interaction of the shock with a particularly dense CSM 
condensation at the far or near end of the shell. The expansion of the SN is 
moderately decelerated (see, Fig [3j lower right panel). 



Bartel 




mas 



1 10 

Time since shock breakout (yr) 



Figure 3. Left: A series of VLBI images of SN 1986J from 1988.7 to 2002.9 
with observing frequency and time after shock breakout (1983.2) indicated. 
The contours are at —8, 8, 11.3, 16, 22.6, 32, ..., and 90.5% of the peak 
brightness. Top right: A 15-GHz image from 2003.5 showing the central com- 
ponent in grey scale and the 5-GHz image from 2002.9 in (slightly different) 
contours, in the same reference frame . Bottom right: The e xpansion of SN 
1986J with a best fit of oc i°- 71±011 (jBietenholz et alll2005h . 



4.3. SN 1979C 



SN 1979C in M100 in Virgo has been studied since 1982, and because of its 
distance of ~15 Mpc could only recentl y be resolved well enough to reveal its 
shell-like structure (Fig. 4, left panel; iBartel fe Bietenhora 120081 ) . The most 
interesting characteristic is that its expansion is consistent with being almost free 
up to t ~20 yr (Fig. 4, right panel). However, significant deceleration may start 
around this time. From fits to the radio lightcurve, a relatively high progenitor 
mass-loss r ate of 1 to 2 x 10~ 4 Mq yr -1 for a w ind velocity of 10 kms -1 was 
estimated ( Weiler et al.lll99ll ; iMontes et al.ll2000l ). However, the kinetic energy 
of the SN explosion would have to be ~ 2 x 10 52 erg to sweep up the mass that 
corresponds to such mass-loss, more than an order of magnitude larger than is 
thought appropriate. The mass- loss rate is therefore likely an order of magnitude 
smaller. Estimates of mass-loss to wind-velocity ratios f rom radio lightcurve 
fittin gs may therefore have to be interpreted with caution (JBartel &: Bietenholz 
20031 ). 
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Figure 4. Left: A 5.0-GHz VLBI image of SN 1979C at t = 25.9 yr. The 
contour levels are at —17, 17, 30, 40, 50 (emphasized), 60, 70, and 90% of the 
peak brightness of 186 /iJy bm _1 . The grey scale, in /uJy bm _1 , is given at 
the top (jBartel fc Bietenholdl20"o1l . Right: The expansion of SN 1979 C with 
a best fit of ex £°- 94±0 - 03 for a shell model (jBartel fc Bietenhold 12001 . 



5. Other optically identified supernovae 

The only other optically identified SN for which detailed images could be ob- 
tained is SN 1987A in the LMC. The images were not obtained with VLBI, 
but rather with the Australia Telescope Compact Array (ATCA). SN 1987A 
is by far the closest SN that has ever been observed with VLBI. Its radio 
lightcurve peaked at ~80 mJy, but when VLBI observations were made in a 
heroic effort t = 5 d after shock breakout, the SN was already completely re- 
solved. Neverthe less, a lower limit o n the expansion velocity of 19,000 km -1 
could be derived (jJauncev et al.lfl988l ). In fact, the SN expanded with a veloc- 
ity of >35,000 kms" 1 and slowed down to 3600 to 5200 kms" 1 for t > 1800 d, 
consistent with the progenitor being a blue supergiant. SN 1987A has a fairly 
round shell or ring- like morphol ogy with a modu lated brightness around the 
ridge ([Manchester et alJl2002l : lGaensler et al.ll2007l ). 

For the other optically identified SNe, no detailed images could be obtained. 
Four of them are Type II. SN 1980K and SN 2004et in NGC 6946 were observed 
only once with VLBI. Th e angular siz e of SN 1980K was determined to be 
<2.0 mas for a shell model ( Bartel|[l985l ) corresponding to an average expansion 
velocity of <11,000 kms^ 1 . For SN 2004et a marginally resolved image, with 
an indication of an asymmetric brightness distribution, and a lower expansion 
velocity limit of 15,700±2000 kms -1 was obtained, and it was inferred tha t 
synchrotron self- absorption is not relevant for this SN ( Marti- Vidal et aUl2007l ). 
For SN 2001g d a size of 0.37±0.08 m as was obtained a couple of years after 
the explosion ( Perez- Torres et al.l 120051 ) . Similarly, for SN 1996cr a preliminary 
size of 10 mas was determined, whic h helped in the int erpretation of this SN as 
possibly being related to SN 1987A (JBauer et al.ll2008l L 
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The remaining four are Type Ib/c SNe . For SN 2007gr, only a detection 
with VLBI was reported (jParagi et al.ll2007l ). For the others, SN 2001em, SN 
2003L, and SN 2008D, upper limits on sizes with corresponding upper limits 
on average expansion velociti es were obtained. Fo r the l atter two SNe, these 
upper limits are relativistic (jSoderberg et al.l 120051 . l2008h . For SN 2003L, at 
92 Mpc the most distant of the SNe in Table 1, the upper limit helped to 
distinguish between absorptio n models, suggesting a (model-dependent) non- 
relativistic nature of the SN ( Soderberg et al]|2005l ). For SN 2001em, a non- 
relativistic limit of the expansion velocity could be determined directly. The 
SN was consid ered to possibly be related to a GRB event but with a misaligned 
relativistic jet ( Granot fe Ramirez-Ruizll2004l ). However, VLBI measurements of 
an expansion velocity of only 5800±10,000 kms -1 an d a proper motion velocity o f 
only 33,000±34,000 kms" 1 disfavor this possibility (|Bietenholz & Bartelll2007h . 



6. Optically unidentified supernovae and supernova remnants 

There are a few dust-obscured starburst galax ies with compact ra dio sources 
that were found to be SN e or SNRs: M82 dBeswick et all 120061 ). Arp 299 
dNeff et all [20041). Arp 220 ([Lonsdale et al.l [20Q6l^). In two other galaxies. NGC 
6240 (iGallimore fc Beswickl 120041 ) . and Mrk 273 (JBondi et al.ll2005h . compact 
radio sources were found with VLBI, which are also thought to be possibly SNe 
or SNRs. The most and longest studied are those in M82. Two sources could 
be imaged at several epochs. The weaker one, 43.31+592, was found to be 
shell-type. It expands with a velocity of 9000-11,000 kms" 1 , which is typical 
for SNe. The strongest of all compact sources in M82, 41.95+575, however, 
is puzzling. The earliest clearly resolved images suggested a sh ell-like struc 



ture, albeit quite elongated with tw o components at opposite ends (jBartel et al. 
19871 : IWilkinson Sz de Bruvnl [19901 ) . Recently, an expansion velocity of 1500- 
2000 kms -1 was determined, which is small for a SN. In addition, the bipolar 
appearance became stronger, and the nature of this source is now less clear 
(JBeswick et al.ll2006l L 



7. Supernova space VLBI in the context of VSOP-2 

From the optically identified SNe studied so far with VLBI, four had flux den- 
sities peaking at ~100 mJy when they became optically thin: SN 1986J, SN 
1987A, SN 1993J, SN 1996cr. Their sizes ranged from 0.55 to 2.8 mas at these 
times (see Table [1]). VLBI observations with VSOP-2 would allow making de- 
tailed images of such SNe starting at the earliest possible time when the optical 
depth decreases below unity and the SN becomes optically thin. "Filming" in 
radio the earliest stages of a SN evolution promises to add invaluable information 
to our understanding of the aftermath of the explosion. 
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